arXiv:1507.07288v2 [nucl-th] 13 Oct 2015 


Nuclear matter equation of state from a quark-model nucleon-nucleon interaction 


K. Fukukawa,*’^ M. Baldo,* G. F. Burgio,* L. Lo Monaco,^ and H.-J. Schulze* 

* INFN Sezione di Catania, Dip. di Fisica, Universita di Catania, Via Santa Sofia 64, F95123 Catania, Italy 
^ Research Center for Nuclear Physics, Osaka University, 10-1 Mihogaoka, Osaka 567-0047, Japan and 
^ Dipartimento di Fisica, Universita di Catania, Via Santa Sofia 64, 1-95123 Catania, Italy 

(Dated: October 14, 2015) 

Starting from a realistic constituent quark model for the nucleon-nucleon interaction, we derive the equation 
of state (EOS) of nuclear matter within the Bethe-Brueckner-Goldstone approach up to three-hole-line level, 
without need to introduce three-nucleon forces. To estimate the uncertainty of the calculations both the gap 
and the continuous choices for the single-particle potential are considered and compared. The resultant EOS 
is compatible with the phenomenological analysis on the saturation point, the incompressibility, the symmetry 
energy at low density and its slope at saturation, together with the high-density pressure extracted from flow 
data on heavy ion collisions. Although the symmetry energy is appreciably larger in the gap choice in the 
high-density region, the maximum neutron star masses derived from the continuous-choice EOS and the gap- 
choice EOS are similar and close to two solar masses, which is again compatible with recent observational data. 
Comparison with other microscopic EOS is presented and discussed. 

PACS numbers: 13.75.Cs, 21.65.Mn, 26.60.Kp 


I. INTRODUCTION 

The nuclear matter equation of state (EOS) is one of the 
central issues in nuclear physics. Its detailed knowledge 
would allow us to connect the data obtained in laboratory ex¬ 
periments on heavy ion collisions (HIC) and the processes that 
characterize the structure and evolution of compact astrophys- 
ical objects like neutron stars (NS) and supernovae. On the 
other hand, laboratory experiments and astrophysical obser¬ 
vations can put meaningful constraints on the nuclear EOS. 
Unfortunately a direct link between phenomenology and the 
EOS is not possible and theoretical frameworks and inputs are 
necessary for the interpretation of the data. In particular the 
EOS above saturation density is much less constrained than 
around or below saturation. 

An intense activity, lasting since several years, has been de¬ 
veloped in analyzing and interpreting the experimental and 
observational data for the purpose of putting severe con¬ 
straints on the EOS and on the corresponding theoretical mod¬ 
els H. It can be recognized from these efforts that a sound 
theoretical and microscopical framework for modeling the 
EOS can be of help for the establishment of firm results on 
the EOS properties. Along these lines a microscopic many- 
body theory based on interactions among nucleons, that stems 
from strong interaction theory, can be of great relevance in re¬ 
ducing the uncertainties that characterize this type of analysis. 

Meson-exchange models of the nucleon interaction have 
been extensively developed since several years, and applied 
to nuclear matter and NS structure within many-body the¬ 
ory. One can mention the variational method 0, the relativis¬ 
tic Dirac-Brueckner-Hartree-Eock (DBHE) 1^, and the non- 
relativistic Bethe-Brueckner-Goldstone (BBG) expansion |0- 
|3l, which have employed different versions of nucleon- 
nucleon interactions inspired by the meson-exchange model. 
In the non-relativistic scheme three-body forces (TBE) have 
been introduced to obtain the correct saturation point of nu¬ 
clear matter iMl. The ambition of all these approaches is to 
devise an elementary interaction among nucleons that is able 


to describe both few-body nuclear systems and nuclear matter 
in agreement with the existing phenomenological data. This 
program has been only partially successful. It turns out in fact 
that it is difficult to reproduce the binding energy of three- and 
four-body systems and at the same time to predict the correct 
saturation point within this scheme. 

More recently the chiral expansion theory to the nucleon 
interaction has been extensively developed Ill2l - I^ . This ap¬ 
proach is based on a deeper level of the strong interaction 
theory, where QCD chiral symmetry is explicitly exploited in 
a low-momentum expansion of the multi-nucleon interaction 
processes. In this approach multi-nucleon interactions arise 
naturally and a hierarchy of the different orders can be es¬ 
tablished. Despite some ambiguity in the parametrization of 
the force and some difficulty in the treatment of many- 
body systems JUt], the method has marked a great progress 
in the microscopic theory of nuclear systems. Indeed it turns 
out that within this class of interactions a compat¬ 

ible treatment of few-nucleon systems and nuclear matter is 
possible. Along the same lines a chiral force has been 
adjusted to reproduce, within a Monte Carlo calculation, the 
binding of ^He and the phase shifts of neutron-alpha scatter¬ 
ing. The same interaction was used to describe neutron mat¬ 
ter. Coupled cluster calculations with chiral forces including 
three-body forces have be en p erformed for neutron and sym¬ 
metric nuclear matter ll^l29l] . and finite nuclei S. 

Another approach inspired by the QCD theory of strong in¬ 
teraction has been developed by a few groups iImI - I^ . In 
this approach the quark degree of freedom is explicitly in¬ 
troduced and the nucleon-nucleon interaction is constructed 
from gluon and meson exchange between quarks, the latter 
being confined inside the nucleons. One of these quark mod¬ 
els of the nucleon-nucleon interaction, named fss2 llj^ [^ . 
is able to reproduce closely the experimental phase shifts and 
the few-body binding energies OSl - l^ . More recently it has 
been shown that the fss2 interaction is able to reproduce 
also the correct nuclear matter saturation point without any 
additional parameter or need to introduce TBE. 
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In this paper we analyze further the fss2 interaction. On 
one hand we compare the results in nuclear matter with ad¬ 
ditional phenomenological constraints, on the other hand we 
extend the EOS based on the fss2 interaction to higher den¬ 
sity and apply it to NS calculations. In this study we use the 
renormalized energy-independent kernel of this model 0. 
Kernels of quark-model nucleon-nucleon interactions are ob¬ 
tained by using the resonating-group method (ROM) for the 
(3^)-(3^) system. Although they are therefore energy depen¬ 
dent, we can eliminate the energy dependence, as we will see 
in the next section. 

This paper is organized as follows. In Section HI] we 
introduce the quark-model nucleon-nucleon interaction fss2 
lIJsl and the formulation for the renormalized RGM ker¬ 
nel. We will briefly review phase shifts and deuteron prop¬ 
erties. The explicit form of the deuteron wave function will 
be given in Appendix In Section |III| we first recapitulate 
the Brueckner-Hartree-Fock (BHF) calculation and the Bethe- 
Faddeev calculation, based on the BBG framework . Then 
the nuclear matter EOS for the fss2 interaction is reported and 
its properties are discussed in relation to the phenomenologi¬ 
cal constraints, also in comparison with some other theoretical 
methods and interactions. Section|IV|is devoted to the calcu¬ 
lations of NS structure. Conclusions are drawn in SectionjVj 


II. QUARK-MODEL BARYON-BARYON INTERACTION 
ESS2 

The fss2 baryon-baryon interaction ll^l44l] is a low-energy 
effective model, which introduces some essential features of 
QCD. The color degree of freedom is explicitly considered 
within the spin-flavor SU (6) approximation and the antisym- 
metrization of quarks is exactly taken into account within the 
framework of RGM. The full model Hamiltonian for the (3^')- 
(3^') system reads 



where m,- and p, are the constituent quark mass and momen¬ 
tum of each particle respectively, and Tq denotes the center- 
of-mass motion. The remaining terms denote the effective 
quark-quark interaction. 

The confinement potential is a phenomenological re¬ 
type potential, which has the favorable feature that it does 
not contribute to the baryon-baryon interactions. We use a 
color analogue of the Fermi-Breit (FB) interaction t/™ with 
explicit quark-mass dependence, motivated by the dominant 
one-gluon exchange process in conjunction with the asymp¬ 
totic freedom of QCD. This includes the short-range repulsion 
and spin-orbit force, both of which are successfully described. 
On the other hand, the medium-range attraction and the long- 
range tensor force, especially those mediated by pions, are 
extremely nonperturbative from the viewpoint of QCD. These 


are therefore most relevantly described by effective meson- 
exc hang e potentials. Compared with the former version FSS 
l!^ 14711 . in which the scalar (S) and the pseudo-scalar (PS) 
nonets are included, the introduction of the vector (V) nonets 
and the momentum-dependent Bryan-Scott term ll4^ greatly 
improves nucleon-nucleon phase shifts li^ and makes fss2 
sufficiently realistic. 

The RGM equation for the relative wave function x(r) is 
given by 

(0 \E - //K{0(3^)0(3^);t (r)}) = 0, (2) 

where (j) {3q) is the three-quark cluster (nucleon) wave func¬ 
tion and is described by a (Os)^ harmonic oscillator with a 
common width parameter. The antisymmetrization operator 
is denoted by si. In our actual calculation, Eq. (|2l) is solved in 
momentum space ii. We rewrite Eq. (|2]i in Schrodinger-like 
form as 


[e-//o-yRGM(e)]z(»’) = 0, (3) 


where e is the two-nucleon energy measured from its thresh¬ 
old in the center-of-mass system and //q is the kinetic energy 
operator. We regard yRGM(£) = iu ~\-G + eK as the non-local 
and energy-dependent potential. Here, Vo is the direct meson- 
exchange kernel, G represents all exchange kernels for the 
kinetic-energy and interaction terms, and K is the exchange 
normalization kernel. 

In the many-body scattering problem, an energy- 
independent potential is desirable, since the energy of a two- 
nucleon pair is not well-defined in the in-medium scattering 
state. The energy dependence of the RGM kernel can be re¬ 
duced by renormalizing the RGM kernel in the following way 
S. We can rewrite Eq. © as 


e _ At-1 /2 (//o + Vo G) ^ 


‘P(r) =0, 


(4) 


where 


N={^{3q)^{3q)\s/\^i3q)^i3q)) (5) 

is the normalization kernel and 'P(r) = xi.'’") is the renor¬ 

malized RGM wave function. If we define the non-local ker¬ 
nel 


W=N-^I^ {Ho + Vd + G)N-^^^-{Ho + Vd + G) (6) 

and the renormalized RGM potential yttGM = 
then Eq. © becomes 


e-z/o-y*^™ 


'P(r) =0. 


(7) 


The detailed procedure to calculate W can be found in Ap¬ 
pendix A of Ref. 0. The asymptotic behavior of ‘P(r) is 
the same as of x(r), because the square root of the normal¬ 
ization kernel approaches unity at large distances. The phase 
shifts derived from ‘P(r) are the same as those from (r). 

In this study, we use a Gaussian representation of the fss2 
potential for numerical simplicity ll^ . This representation 
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FIG. 1. (Color online) Phase shifts of the energy-independent Gaussian fss2 potential m to the energy 7[ab < 350 MeV (black solid curves), 
compared with the results for the Argonne VI8 1^ (red dashed curves) and CD-Bonn 15^ (blue short-dashed curves) potentials. The symbols 
show results from the Nijmegen multi-energy phase shift analysis (PSA) 03,1551] and Arndt’s tip phase shift analyses l5^ - [^ . 


conserves the non-local feature of the potential accurately and 
reproduces the original phase shifts within an accuracy bet¬ 
ter than 0.1° for almost all energies and partial waves. The 
energy dependence of the Gaussian-represented potential is 


eliminated by the above-mentioned method. The phase shifts 
for the energy-independent version of fss 2 are shown in Fig.[T] 
One can hnd some dehciencies in the ^f ’2 ^^2 3 phase 

shifts. As already pointed in Ref. this is probably re- 
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TABLE I. Neutron-proton phase shifts and proton-proton *So phase shifts (in degrees), and values of the fss2 model with respect to the 
Nijmegen PSA. 


Tiab (MeV) 

0.383 

1 

5 

10 

25 

50 

100 

150 

215 

320 

No. of data 

144 

68 

103 

290 

352 

572 

399 

676 

756 

954 

AX^ 

3034 

2100 

274 

710 

919 

6804 

1409 

3913 

1464 

2637 

'Sq{pp) 

14.519 

33.097 

55.765 

56.601 

50.693 

41.324 

27.350 

16.822 

6.143 

-7.090 

^Soinp) 

54.592 

62.220 

64.223 

60.862 

52.234 

41.916 

27.496 

16.899 

6.130 

-7.105 

^Si 

159.429 

147.838 

118.380 

102.881 

80.956 

63.014 

43.191 

30.495 

18.540 

4.737 

£\ 

0.027 

0.103 

0.653 

1.116 

1.678 

1.915 

2.209 

2.684 

3.553 

5.344 

^Pq 

0.048 

0.195 

1.747 

3.885 

8.546 

11.250 

9.038 

4.018 

-3.147 

-14.100 

^P\ 

-0.030 

- 0.122 

-1.033 

-2.249 

-5.228 

-8.685 

-13.448 

-17.273 

-21.786 

-28.502 


-0.049 

-0.193 

-1.545 

-3.183 

-6.703 

-10.267 

-14.826 

-18.392 

-22.502 

-28.421 

^Pl 

0.005 

0.020 

0.237 

0.682 

2.573 

6.252 

12.412 

15.901 

17.485 

16.592 

£2 

- 0.000 

- 0.002 

-0.057 

-0.207 

-0.824 

-1.767 

-2.854 

-3.218 

-3.211 

-2.865 

^Di 

- 0.001 

-0.005 

-0.183 

-0.679 

-2.821 

-6.519 

-12.427 

-16.596 

-20.160 

-22.872 

3^2 

0.001 

0.006 

0.221 

0.842 

3.671 

8.825 

17.091 

22.260 

22.562 

26.213 

IO 2 

0.000 

0.001 

0.045 

0.162 

0.635 

1.473 

3.294 

5.304 

7.816 

10.677 

3^3 

0.000 

0.000 

0.001 

0.001 

0.008 

0.192 

1.287 

2.975 

5.361 

8.425 

£3 

0.000 

0.000 

0.013 

0.081 

0.554 

1.625 

3.570 

5.008 

6.250 

7.271 


0.000 

0.000 

0.002 

0.013 

0.101 

0.319 

0.718 

0.978 

1.076 

0.595 

3F3 

- 0.000 

- 0.000 

-0.005 

-0.033 

-0.230 

-0.703 

-1.612 

-2.315 

-2.940 

-3.471 

'F 3 

- 0.000 

- 0.000 

- 0.011 

-0.066 

-0.424 

-1.156 

-2.344 

-3.194 

-4.028 

-5.094 


- 0.000 

0.000 

0.000 

0.001 

0.017 

0.087 

0.366 

0.816 

1.642 

3.455 

£4 

- 0.000 

- 0.000 

- 0.000 

-0.004 

-0.047 

-0.191 

-0.540 

- 0.866 

-1.223 

-1.642 

3G3 

- 0.000 

- 0.000 

- 0.000 

-0.004 

-0.054 

-0.264 

-0.982 

-1.882 

-3.109 

-4.889 

3G4 

0.000 

0.000 

0.001 

0.014 

0.171 

0.726 

2.191 

3.673 

5.423 

7.716 

‘G4 

0.000 

0.000 

0.000 

0.003 

0.038 

0.147 

0.386 

0.606 

0.886 

1.413 


0.000 

- 0.000 

- 0.000 

- 0.000 

-0.009 

-0.052 

-0.189 

-0.316 

-0.399 

-0.279 

£5 

0.000 

0.000 

0.000 

0.002 

0.037 

0.206 

0.733 

1.300 

1.983 

2.917 


TABLE 11. Comparison of the deuteron properties with the predictions using fss2, Bonn C, CD-Bonn, and experimental data. 



fss2-Gauss 

fss2 (isospin) 1381 

Bonn C 

CD-Bonn 

Expt. 

£d (MeV) 

2.2206 

2.2250 

fitted 

fitted 

2.224644±0.000046 

Rrms (fm) 

1.961 

1.960 

1.968 

1.966 

1.971 ±0.006 162-641 

Qd (fm^) 

0.270 

0.270 

0.281 

0.270 

0.2859±0.0003 165. 66 ] 

77 = Ad/A s 

0.0252 

0.0253 

0.0266 

0.0256 

0.0256±0.0004 [W] 

Pd (%) 

5.52 

5.49 

5.60 

4.85 



lated to the problem of the balance of central and LS force in is taken into account by a reduction factor for the coupling 


the short-range region. Appreciable deviations appear also in 
the ^Di and channels in the higher-energy region. This 
implies that improvements of the tensor force are desirable in 
future refinements of the interaction. 

Examining with respect to the PSA is a good test to 
see how strict the interaction describes the observables. In 
this study, we worked in the isospin basis and with the cut-off 
Coulomb force 


Vc{k) 


1 - cos(k/?cut) 

I? 


( 8 ) 


and extracted the nuclear phase shifts by using the Vincent- 
Phatak method llb^ . This method give s stable phase shifts 
with respect to the change of Rcut ll^ . which we take as 
Rcut = 10 fm. 

Another important factor in considering the pp phase shifts 
is the charge-independence breaking (CIB). The CIB effect 


constant of the scalar-singlet meson, which is determined to 
minimize x'^ ll^ Idd]] . Our value of the reduction factor is 
0.9932, which is quite close to 0.9934 used in Ref. Idd]] . We 
report in Tab. Uthe neutron-proton (iip) and the proton- 
proton (pp) phase shifts with their Ax'^ values with respect 
to the Nijmegen PSA The Gaussian fss2 potential 

gives X^/Ndata = 6.34. Although the phase shifts are overall 
well reproduced, as we saw in Fig.[T] we have larger x^ values 
mainly due to the ^P 2 , ^D 2 , and partial waves in the high- 
energy region. Moreover, the constraint for the pp phase 
shifts in the low-energy region is very severe. Those phase 
shifts given by the Nijmegen PSA are 14.609° and 32.688° 
at Ep = 0.383 MeV and 1 MeV, respectively if^ . A more 
developed treatment is desirable, because the difference from 
the Nijmegen PSA should be much less than 0.1°. 

When the deuteron properties were determined in Ap- 
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FIG. 2. (Color online) The deuteron wave functions u{r) = Moa(^) 
and w{r) = U 2 a{r) predicted by the energy-independent Gaussian po¬ 
tential based on fss2 (black solid curve), compared with CD-Bonn 
(blue dashed) j53|]. 



FIG. 3. (Color online) As Fig.|2l but for the deuteron wave functions 
in momentum space qu{q) = qfoaiq) and = qflaiq)- 


pendix B of Ref. the authors first solved Eq. (O for 

X(r) and then obtained the renormalized relative wave func¬ 
tion ‘E(r) = %{'<')■ For the renormalized potential, we 

can directly solve Eq. dUi for 'E(r). Because the procedure to 
obtain the potential and the deuteron wave function is differ¬ 
ent, it is necessary to reexamine the deuteron properties. The 
detailed prescription is explained in Appendix A. 

Eigures |2] and [3 show the deuteron wave functions by the 
energy-independent fss2 potential in coordinate and momen¬ 
tum space, respectively. The wave functions are indistinguish¬ 
able from those of Ref. . It can also be seen from Table 
im that the renormalization of the kernel does not change the 


(«) 

o—o 


(l>) 



FIG. 4. Different Goldstone diagrams contributing to the nuclear 
matter EOS. Diagrams (a) and (b) correspond to the BHF calculation. 
The sum of the other diagrams gives the three-hole-line contribution. 
For more details, please see the text. 


deuteron properties. Our calculation slightly underpredicts 
the quadrupole moment, as other potentials do. According 
to Refs. IITllIT^ . the correction due to the meson-exchange 
current is typically about 0.01 fm^. However, even taking into 
account this correction, there is still a small discrepancy. 


III. NUCLEAR MATTER EOS WITHIN THE 
BETHE-BRUECKNER-GOLDSTONE APPROACH 

A. Sketch of the approach 

The basis of the BHF calculation is the Bethe-Goldstone 
equation for the G matrix, 

(12|G(a))|34)^ = (12|y|34)^ 

+ ^(12|y|56)—-(56|G(®)|34)^ , (9) 

te (O-es-ee 

where the multi-indices 1,2 ,... include the momentum and 
the spin-isospin variables of the particles, 112)^ = 112) — |21), 
y is the bare nucleon-nucleon interaction, (O is the starting 
energy, and = 0(^5 — — is the two-particle 

Pauli operator, where k^p are the Fermi momenta of the nu¬ 
cleons. In this study, we use the angle-averaged form of 
the Pauli operator and the two-particle intermediate energy 
es + e^, in order to avoid the complex coupling of the angular- 
momentum quantum numbers 11731] . 

The single-particle (s.p.) energy is defined by e, = 
kf /2M]^ + Ui{ki), where is the nucleon mass. The aux¬ 
iliary s.p. potentials G, (k,) are self-consistently determined by 
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TABLE III. Energies per nucleon (in MeV) of the various THL contributions to the SNM EOS for different Eermi momenta kp in fm^*. The 
continuous choice is adopted. The baryon density p = 2kp/{37t^) in fm^^ is also shown. T + E 2 corresponds to the BHF calculations, B 
is the bubble diagram, BU is the (/-insertion diagram, R is the ring diagram, H indicates the higher-order diagrams, and £3 is the total THL 
contribution. 


kp 

P 

T + E 2 

B 

BU 

R 

H X 10^ 

£3 

T + £2 + £3 

1.1 

0.090 

-17.07 

-7.14 

11.02 

-0.77 

186.8 

3.30 

-13.77 

1.2 

0.117 

-19.58 

-6.17 

11.46 

-1.27 

159.3 

4.22 

-15.36 

1.3 

0.148 

-21.98 

-4.37 

11.68 

-1.63 

113.8 

5.79 

-16.19 

1.4 

0.185 

-24.21 

-2.33 

12.48 

-1.87 

77.22 

8.35 

-15.86 

1.5 

0.228 

-26.91 

0.80 

12.93 

-1.99 

48.35 

11.79 

-14.30 

1.6 

0.277 

-27.41 

4.93 

13.91 

-2.13 

26.58 

16.73 

-10.67 

1.7 

0.332 

-27.90 

10.76 

14.83 

-2.05 

10.88 

23.56 

-4.34 

1.8 

0.394 

-27.33 

18.91 

15.61 

-1.93 

1.19 

32.59 

5.27 

1.9 

0.463 

-25.57 

29.85 

16.21 

-1.75 

-3.75 

44.31 

18.74 

2.0 

0.540 

-22.25 

43.71 

15.09 

-1.42 

2.49 

57.39 

35.14 

2.1 

0.626 

-17.08 

61.42 

13.52 

-1.29 

0.80 

73.66 

56.57 

2.2 

0.719 

-9.73 

83.58 

10.47 

-1.19 

- 0.11 

92.86 

83.13 

2.3 

0.822 

- 0.01 

111.12 

6.15 

- 1.11 

- 0.21 

116.16 

116.16 


TABLE IV. As Tablelllll but with the gap choice for the s.p. potential. 

kp 

p 

T+E 2 

B 

R 

Hx 10^ 

£3 

£ + £2 + £3 

1.1 

0.090 

-11.54 

-0.48 

-1.00 

65.07 

-1.42 

-12.95 

1.2 

0.117 

-13.45 

0.08 

-1.09 

49.27 

-0.96 

-14.41 

1.3 

0.148 

-15.26 

1.02 

-1.18 

33.16 

-0.13 

-15.39 

1.4 

0.185 

-16.86 

2.48 

-1.22 

20.72 

1.28 

-15.59 

1.5 

0.228 

-18.12 

4.64 

-1.25 

11.66 

3.40 

-14.72 

1.6 

0.277 

-18.88 

7.82 

-1.29 

5.45 

6.54 

-12.34 

1.7 

0.332 

-19.00 

12.40 

-1.26 

0.99 

11.14 

-7.85 

1.8 

0.394 

-18.31 

18.96 

-1.24 

-1.99 

17.72 

-0.59 

1.9 

0.463 

-16.65 

28.03 

-1.19 

-3.96 

26.84 

10.19 

2.0 

0.540 

-13.90 

39.71 

-1.09 

-1.80 

38.62 

24.73 

2.1 

0.626 

-9.90 

55.18 

-1.06 

-1.83 

54.12 

44.22 

2.2 

0.719 

-4.53 

75.75 

-1.95 

-1.76 

74.71 

70.17 

2.3 

0.822 

2.29 

102.63 

-1.06 

-1.56 

101.57 

103.85 


the on-shell G matrix elements, with Eq. (|9]l, 


asymmetric nuclear matter by the parabolic approximation. 


Ui{ki) = Y. {iJ\G{ei + ej)\ij)^. (10) 

J<kF 


The s.p. potential U{k) can be chosen in various ways. Our 
investigations will be carried out for two somewhat oppo¬ 
site choices: the continuous choice and the gap (or standard) 
choice. In the continuous choice, Eq. (fTOl) is solved for all 
k, while U{k) =0 is assumed for ^ in the gap choice. 
The detailed procedure for BHE calculations is presented in 
Refs. El and S. 

In this study, we calculate the binding energy per particle 
for symmetric nuclear matter (SNM) £’/A(p,Xp = 0.5) and 
pure neutron matter (PNM) ElA{p,Xp = 0), where Xp = Pp/p 
is the proton fraction. In these cases the energy per particle 
from the two-hole-line contribution is given by 


E\ 3 k} 
a) 2 5 2Mn 


1 


L u{k). 

k<kf 


( 11 ) 


At a given baryonic density p, we will approximate E jA for 


l[p^xp) = {\-fi)j 


ip) + P- 


SNM 


(P)> 

PNM 


( 12 ) 


where j3 = (1 — 2xp)^, which has been verified to be a good 
approximation within the BHE approach ItI . 

The two-hole-line and three-hole-line (THL) diagrams are 
depicted in Eig. |4l where the wavy line denotes the G ma¬ 
trix. Eigures Ufa) and |4|b) are the above-mentioned BHE 
direct (Hartree) and exchange (Lock) diagram, respectively. 
As for the THL calculations, we closely follow the method 
described in detail in Ref. H. The full scattering process 
of three particles that are virtually excited above the Eermi 
sphere can be calculated by solving the Bethe-Eaddeev equa¬ 
tion for the in-medium three-body scattering matrix as 
depicted in Eig. |4tf) |0, Esl I3. Eor computational conve¬ 
nience, the lowest-order contribution in the G matrix, namely 
diagram Etc), is calculated separately. EigureEfc) is known 
as the “bubble” diagram, and Eig. E^d) is the corresponding 
(/-insertion diagram. Note that (/-insertion diagram vanishes 
in the gap choice, because U{k) =0 is assumed for k > kp. 
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TABLE V. As TableHni but for PNM. 


kp 

P 

T + E2 

B 

BU 

R 

H X 10^ 

£3 

T + E 2 +E 2 

1.0 

0.034 

5.18 

-0.84 

1.17 

-0.01 

2.63 

0.33 

5.50 

1.1 

0.045 

5.96 

-0.85 

1.38 

-0.11 

-0.04 

0.42 

6.38 

1.2 

0.058 

6.83 

-0.76 

1.53 

-0.15 

-3.74 

0.62 

7.45 

1.3 

0.074 

7.66 

-0.63 

1.84 

-0.15 

-7.02 

1.04 

8.71 

1.4 

0.093 

8.55 

-0.38 

2.14 

-0.15 

-9.34 

1.60 

10.15 

1.5 

0.114 

9.56 

0.07 

2.41 

-0.13 

-10.59 

2.34 

11.90 

1.6 

0.138 

10.71 

0.74 

2.81 

-0.10 

-11.32 

3.44 

14.14 

1.7 

0.166 

12.10 

1.79 

3.25 

-0.09 

-11.17 

4.94 

17.04 

1.8 

0.197 

13.87 

3.40 

3.66 

-0.09 

-10.39 

6.96 

20.83 

1.9 

0.232 

16.05 

5.69 

4.11 

-0.09 

-8.81 

9.70 

25.75 

2.0 

0.270 

18.76 

8.74 

4.30 

-0.05 

-5.54 

12.98 

31.74 

2.1 

0.313 

22.23 

12.81 

4.39 

-0.07 

-4.17 

17.13 

39.36 

2.2 

0.360 

26.61 

18.92 

4.14 

-0.08 

-2.99 

22.97 

49.59 

2.3 

0.411 

32.02 

24.66 

3.66 

-0.10 

-2.09 

28.22 

60.23 

2.4 

0.467 

38.53 

33.09 

2.76 

-0.10 

-1.33 

35.74 

74.28 

2.5 

0.528 

46.42 

43.30 

1.18 

-0.11 

-0.75 

44.36 

90.79 

2.6 

0.594 

55.86 

55.46 

-1.08 

-0.13 

-0.39 

54.25 

110.11 

2.7 

0.665 

66.93 

70.50 

-4.14 

-0.14 

-0.15 

66.22 

133.15 

2.8 

0.741 

79.78 

89.54 

-8.23 

-0.15 

-0.09 

81.16 

160.94 

2.9 

0.824 

94.52 

110.30 

-13.50 

-0.27 

0.42 

96.53 

191.05 


TABLE VI. As Table Ivl but with the gap choice for the s.p. potential. 

kf 

P 

T+E 2 

B 

R 

H X 10^ 

£3 

T + £2 + £3 

1.0 

0.034 

5.79 

-0.23 

-0.05 

1.03 

-0.28 

5.51 

1.1 

0.045 

6.70 

-0.18 

-0.08 

-0.82 

-0.26 

6.43 

1.2 

0.058 

7.63 

-0.09 

-0.10 

-2.68 

-0.18 

6.51 

1.3 

0.074 

8.62 

0.07 

-0.09 

-4.21 

-0.03 

8.59 

1.4 

0.093 

9.69 

0.32 

-0.09 

-5.10 

0.23 

9.91 

1.5 

0.114 

10.88 

0.71 

-0.07 

-5.42 

0.63 

11.50 

1.6 

0.138 

12.25 

1.30 

-0.06 

-5.56 

1.24 

13.49 

1.7 

0.166 

13.88 

2.19 

-0.05 

-5.45 

2.14 

16.02 

1.8 

0.197 

15.85 

3.53 

-0.05 

-5.16 

3.47 

19.32 

1.9 

0.232 

18.24 

5.43 

-0.06 

-4.60 

5.37 

23.61 

2.0 

0.270 

21.15 

7.95 

-0.05 

-3.11 

7.90 

29.06 

2.1 

0.313 

24.69 

11.34 

-0.06 

-2.41 

11.28 

35.97 

2.2 

0.360 

28.96 

15.89 

-0.08 

-1.79 

15.19 

44.77 

2.3 

0.411 

34.40 

21.76 

-0.09 

-1.26 

21.67 

55.71 

2.4 

0.467 

40.02 

29.50 

-0.10 

-0.83 

29.40 

69.43 

2.5 

0.528 

47.00 

39.64 

-0.12 

-0.47 

39.53 

86.52 

2.6 

0.594 

55.02 

52.63 

-0.14 

-0.24 

52.49 

107.51 

2.7 

0.665 

64.13 

69.74 

-0.16 

-0.06 

69.58 

133.71 

2.8 

0.741 

74.37 

93.04 

-0.19 

-0.05 

92.85 

167.22 

2.9 

0.824 

85.70 

123.18 

-0.21 

0.08 

122.97 

208.67 


Figure |4te) is the “ring” diagram, which is responsible for 
long-range correlations in nuclear matter. An indication of 
convergence of the BBG expansion is the possible small size 
of the THL contribution with respect to the two-hole-line con¬ 
tribution. 


B. Numerical results 

We report in Tables HIIllVIl the contributions of each dia¬ 
gram to the EOS for SNM and PNM in the continuous and 


gap choice, respectively. The slightly different results with 
respect to Ref. are due to the more refined momentum 
grid we used in the present calculations. This was relevant 
at higher density in order to obtain convergence in the BHF 
iteration procedure. We divide the space of relative momen¬ 
tum q into the two domains, [0,a] and [a,°o] and apply the 
Gauss-Legendre quadrature to each part in solving the Bethe- 
Goldstone equation. The mapping q = a + tan[(1 +x)n/A], 
where x are the nodes of the Gauss-Legendre quadrature, is 
used for the second part as in Ref. lld^ . In this study, a = 6 fm 
is adopted and we take 70 points in the first section and 30 
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TABLE VII. The fitted coefficients in Eq. ( fTst for SNM and PNM 
and continuous (C) or gap (G) choice, together with those obtained 
for other interactions. 



a 

b 

c 

d 

SNM,C 

-515.4 

692.6 

1.30 

2.4 

SNM,G 

-113.6 

296.4 

1.91 

-6.2 

PNM,C 

-103.5 

355.3 

1.48 

8.3 

PNM,G 

60.2 

252.3 

2.50 

3.4 

SNM,APR 

-101.5 

333.9 

2.14 

-4.8 

SNM,DBHP 

-422.8 

711.3 

1.56 

7.3 

SNM,V18-l-TBP 

-123.2 

407.9 

2.38 

0 . 

SNM,V18-l-UIX 

-137.0 

308.0 

1.82 

-5.0 

PNM,APR 

76.1 

256.5 

2.71 

3.6 

PNM,DBHP 

-230.3 

715.5 

1.58 

9.5 

PNM,V18-l-TBF 

55.9 

532.3 

2.68 

0 . 

PNM,V18-l-UIX 

11.0 

309.0 

1.95 

6.0 


points in the second section. All nucleon-nucleon channels up 
to the total angular momentum 7 = 8 were considered. In each 
iteration, the s.p. potential was calculated self-consistently up 
to A:max = 7.5 fm^^ with a grid step of 0.1 fm^^ After 30 
iterations, a convergence within few keV was reached in all 
calculations at kp <3.0 fm^*. 

After adding the THL contributions, we fitted the calculated 
EOS, both for SNM and PNM, by an analytic form with four 
parameters 

^{p)=ap+bp‘' + d. (13) 

The values of the fitted parameters are listed in Table VII. 
These fitted EOS are valid for not too low density p > 
0.1 fm^^ and must not be extrapolated to zero density. In 
the continuous choice around saturation the analytic form is 
very close to the calculated points, and one can extract the 
saturation energy {E/A)o = —16.3 MeV at po = 0.157 fm^^ 
and an incompressibility K = 219 MeV, while with the gap 
choice one finds {E/A)q = —15.6 MeV, po = 0.170 fm^^, 
and /T = 185 MeV. This indicates the uncertainty of the EOS 
around saturation. At higher density the fit is less precise, but 
the deviation does not exceed 1.2 MeV even at the highest 
density in the continuous choice, which is more than enough 
for NS calculations. In the gap choice, it is not so easy to de¬ 
scribe the EOS by one single analytic expression from very 
low to high density as in the continuous choice, but the re¬ 
ported values of the parameters give a good fit at high density, 
which is useful for the NS study. 

The fss2 EOS for SNM and PNM are reported in Pig.|5] in 
comparison with the corresponding EOS from some other ap¬ 
proaches. The latter have been selected from the ones that are 
able to reproduce the saturation point within the phenomeno¬ 
logical uncertainty. The comparison of the different EOS for 
SNM shows a substantial agreement up to about 0.5 fm^^, 
while at higher density both the variational calculation (AP^ 
of Ref. 101 and the relativistic DBHE calculation of Ref. ||3D 
indicate a stiffer trend. As for non-relativistic BHE calcu¬ 
lations, also the EOS obtained with the “microscopic” TBE 
of Refs, in [HI] is appreciably stiffer in that density region. 


while the BHE with the Urbana model jHHHj for the TBE pro¬ 
duces an EOS in substantial agreement with the fss2 EOS. A 
similar trend is present in PNM. 

In Pig. |5] the EOS is reported in the density range relevant 
for NS calculations. Prom Tables HIIllVIl one can notice that 
at the highest densities the THL contribution can be larger 
than the two-hole-line one. This can make the convergence of 
the BBG expansion at least questionable. However, it turns 
out that the two-hole-line interaction part is quite limited be¬ 
cause of the strong compensation between negative and posi¬ 
tive channel contributions, and actually it is decreasing in the 
highest density region. On the contrary, the dominant contri¬ 
bution to the THL interaction term is coming from the repul¬ 
sive “bubble” diagram. This consideration can suggest that 
the convergence might still be present, but of course it cannot 
give a strong argument in support of it. 

However, in addition to that, one can see from Tables m 
|yl]that the three-body scattering processes, described by the 
scattering matrix (column “//”), give a negligible contri¬ 
bution. It can then be expected that the four-body scattering 
processes will also be negligibly small. The fourth-order di¬ 
agrams apart from the four-body scattering processes, have 
been estimated in Ref. up to about three times saturation 
density for the Reid soft-core NN interaction iFtsI] and found 
to be quite small. We will assume that this is also true in our 
case, even for higher density, at least approximately enough 
for NS calculations. 

As a check of this assumption we will confront the EOS 
with known phenomenological constraints. The higher den¬ 
sity part of the EOS, needed for NS calculations, can be seen 
as an extrapolation from the lower one, which can be validated 
from the comparison with astrophysical observations and lab¬ 
oratory experiments on heavy-ion collisions. It is clear that 
the main theoretical uncertainty on the EOS is indeed coming 
from the many-body calculations at higher density. Unfortu¬ 
nately it is difficult to get a quantitative estimate of this uncer¬ 
tainty without a firm limit on the higher-order contributions 
beyond the THL. 


C. Comparison with phenomenology 

We will now confront the EOS with a set of phenomeno¬ 
logical constraints in order to assess its reliability. Possible 
tests of the EOS have been devised from experiments on HIC, 
which have been performed in the last two decades at ener¬ 
gies ranging from few tens to several hundreds MeV per nu¬ 
cleon. It can be expected in fact that in HIC at large enough 
energy nuclear matter is compressed and that the two part¬ 
ners of the collisions produce flows of matter. In principle the 
dynamics of the collisions should be connected with the prop¬ 
erties of the nuclear medium EOS and its viscosity. In the so- 
called “multifragmentation” regime, after the collision numer¬ 
ous nucleons and fragments of different sizes are emitted, and 
the transverse flow, which is strongly affected by the matter 
compression during the collision, can be measured. Based on 
numerical simulations, it was proposed in Ref. that any 
reasonable EOS for SNM should pass through a phenomeno- 
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FIG. 5. (Color online) The energy per particle of SNM (upper panel) 
and PNM (lower panel) for several EOS as a function of baryon den¬ 
sity. The curves FSS2 (CC) and FSS2 (GC) show the calculated 
fss2 EOS in continuous and gap choice, respectively. The curve 
Avig-l-UIX refers to the BHE calculation using Argonne V18 plus 
Urbana IX TBF from Ref. |13,Sl> the curve APR is the variational 
calculation from Re f. Ip . the curve Avig-l-micro TBF is the BHF cal¬ 
culation from Ref. fTo niTIl . and the curve DBHF is the relativistic 
DBHF calculation from Ref. d]. 



FIG. 6. (Color online) The pressure of symmetric matter for several 
EOS. The larger (yellow) and smaller (violet) bands represent the 
phenomenological constraints from experimental data. See text for 
details. 



FIG. 7. (Color online) The symmetry energy for several EOS. 


logical region in the pressure vs. density plane. 

The plot is reproduced in Fig. |6] where a comparison with 
the same microscopic calculations is made. The larger (yel¬ 
low) dashed box represents the results of the numerical simu¬ 
lations of the experimental data discussed in Ref. if?^ . and the 
smaller (violet) one represents the constraints from the exper¬ 
imental data on kaon production 1^ . The fss2 EOS is in any 
case fully compatible with the phenomenological constraints. 
This is true also for the other selected EOS, with the possi¬ 
ble exception of the DBHE one, which appears too repulsive 
at higher density. The analysis indicates that the EOS at low 
density must be relatively soft. 

Another gross property of the nuclear EOS, which plays a 
decisive role in NS calculations, is the symmetry energy as a 
function of density S{p), especially at the high density typical 
of the NS inner core. It can be expressed in terms of the energy 
per particle between PNM and SNM, 

S{P) = -i^^(p,0.5) ^ |(p,0) - |(p,0.5), (14) 

and is reported in Eig. |2]for the considered set of EOS. A 
large spread of values is present for densities above saturation. 
In comparison with the other EOS, the fss2 EOS appears in 
the region of lower values (“iso-soft” EOS), but the gap and 
continuous choices show an appreciable discrepancy at higher 
density. As we will see, the stiffness of NS matter shows a 
reduced spread of values, since the beta-equilibrium condition 
produces a compensation effect due to the interplay of the size 
of the symmetry energy and the stiffness of the EOS. 

The symmetry energy up to saturation density has been con¬ 
strained in Ref. M by analyzing the data on isobaric analog 
states as well as on the neutron skin in a set of nuclei. In Eig. [8] 
the larger (yellow) band indicates the constraint coming from 
the analog states, while the more restricted region bounded by 
the full (red) line is obtained if also the neutron skin data are 
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FIG. 8. (Color online) The symmetry e nerg y at low density. The 
yellow band represents recent constraints 18111 and the red line shows 
the region restricted by the neutron skin data, whereas the different 
curves are the results of the microscopic many-body methods. 


added. The fss2 EOS is consistent throughout the constrained 
regions. The other EOS look also consistent with the con¬ 
straints, with the possible exception of the BHE calculation 
with microscopic TBE 

Another parameter which characterizes the symmetry en¬ 
ergy is its slope at saturation, usually embodied in the quan¬ 
tity L = 3p((95'/(9p)|p=p(,. In Eig. |9]this parameter is dis¬ 
played versus the value of the symmetry energy at saturation 
5o = S{po), which has been widely discussed in Ref. lls^l . The 
different boxes indicate several constraint regions obtained in 
different phenomenological analyses. The dashed (blue) band 
represents the constraint coming from experimental data on 
HIC, obtained from the neutron and proton spectra from cen¬ 
tral collisions for *^^SnH-*^^Sn and *^^SnH-**^Sn reactions at 
50 MeV/A iHl]. At the same incident energy, isospin diffu¬ 
sion was investigated. We remind that isospin diffusion in 
HIC depends on the different N/Z asymmetries of the in¬ 
volved projectiles and targets, hence it is used to probe the 
symmetry energy ifsil - f^ . The full black circle shows the 
results from isospin diffusion observables measured for col¬ 
lisions at a lower beam energy of 35 MeV per nucleon 0. 
Transverse collective flows of hydrogen and helium isotopes 
as well as intermediate mass fragments with Z <9 have also 
been measured at incident energy of 35 MeV/A in ^*'ZnH-^‘'Zn, 
^^Zn+^^Zn, ®^NiH-®^Ni reactions and compared to transport 
calculations. The analysis yielded values denoted by the full 
black squares . 

The box labeled by ERDM (finite-range droplet model) rep¬ 
resents a refinement of the droplet model 11 ^ . and includes 
microscopic “shell” effects and the extra binding associated 
with N — Z nuclei. The ERDM reproduces nuclear binding 
energies of known nuclei within 0.1%, and allows determina¬ 
tion of both So = 32.5 ± 0.5 MeV and L = 70 ± 15 MeV. 

The other boxes in Eig. |9] represent experimental data ob- 



FIG. 9. (Color online) L = 3p{dS/dp)\p =p„ versus symmetry en¬ 
ergy at saturation 5o, predicted by several EOS. See text for details 
of the experimental constraints. 


tained from measurements of the neutron skin thickness. In 
light nuclei with N k,Z, the neutrons and protons have similar 
density distributions. With increasing neutron number V, the 
radius of the neutron density distribution becomes larger than 
that of the protons, reflecting the pressure of the symmetry en¬ 
ergy. The measurement of the neutron skin thickness is made 
on the stable nucleus which has a closed neutron shell 

with V = 126 and a closed proton shell with Z = 82, hence 
it is very asymmetric and the neutron skin is very thick. The 
possibility of measurements of the neutron radius in 
the experiment PREX at Jefferson Laboratory has been widely 
discussed ll^ . The experiment should extract the value of the 
neutron radius in from parity-violating electron scatter¬ 
ing. However, the experimental signature is very small, and 
the extracted thickness has a large statistical uncertainty. In 
the next few years, a second experimental run for PREX could 
reduce this large uncertainty 11 ^ . 

Recent experimental data obtained by Zenihiro et al. ll93] 
on the neutron skin thickness of 208p|^ deduced a value of 
= 0.21 l)')QQg 3 fm. Erom the experiments constraints 
on the symmetry energy were derived, and these are plot¬ 
ted in Eig. |9]as the short-dashed blue rectangular box labeled 

Pb(p,p'). 

Last, we mention the experimental data on the Pygmy 
Dipole Resonance (PDR) in very neutron-rich nuclei such as 
®^Ni and ^^^Sn, which peaks at excitation energies well below 
the Giant Dipole Resonance, and exhausts about 5% of the 
energy-weighted sum rule roll . In many models it has been 
found that this percentage is linearly dependent on the slope 
L of the symmetry energy. Values of L = 64.8 ± 15.7 MeV 
and = 32.2 ±1.3 MeV were extracted in Ref. lodl] . using 
various models which connect L with the neutron skin thick¬ 
ness. Those constraints are shown in Pig.|9]as a long-dashed 
rectangle with the label PDR. 

It is not clear to what extent all these constraints are com- 
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FIG. 10. (Color online) Proton fraction Xp (upper panel), pressure 
P (middle panel), and speed of sound in units of c (lower panel) of 
beta-stable nuclear matter for several EOS as a function of baryon 
density p. 


patible among each other, but it looks that most of the EOS 
provide values consistent with the general trend, including the 
fss2 EOS. 


IV. NEUTRON STAR STRUCTURE 


In order to study the structure of neutron stars, we have to 
calculate the composition and the EOS of cold, neutrino-free, 
catalyzed matter. We require that the neutron star contains 
charge-neutral matter consisting of neutrons, protons, and lep¬ 
tons (e^, p in beta equilibrium, and compute the EOS in the 
following standard way ll^ 1^ : The Brueckner calculation 
yields the energy density of b ary on/lepton matter as a function 
of the different partial densities. 


^{Pnj Ppi Pe^ Pjl) — {PiP^n T Ppf^p) T (Pn T Pp) ^ (Pn j Pp) 


+ Ptimp 


1 

2m a 


471^ 


(15) 


where we have used ultrarelativistic and nonrelativistic ap¬ 
proximations for the energy densities of electrons and muons. 


respectively. In this study, we adopted the parabolic ap¬ 
proximation for E/A, Eq. (fT^ . Knowing the energy den¬ 
sity Eq. dni), the various chemical potentials (of the species 
i = n,p,e,ii) can be computed straightforwardly, 

de 

P/ = 3— , ( 16 ) 

dpi 

and the equations for beta equilibrium. 

Pi — ^iPn iJiPe ^ ( 17 ) 


{bi and p, denoting baryon number and charge of species i) 
and charge neutrality, 

Y.P^qi = 0, (18) 


allow one to determine the equilibrium composition {p,} at 
given baryon density p and finally the EOS, 


P{p)=p 


2 d e({p,(p)}) de 


dp 


= p—-e = pp„-e. 
dp 


(19) 


In Eig. [T0]we report the proton fraction Xp, the pressure P, 
and the sound velocity. 


Cs = 



( 20 ) 


as a function of total baryon density in NS matter. The sound 
velocity can be used as a further test of a given EOS, since it 
should not exceed the speed of light c. One can see that the 
fss2 EOS is becoming superluminal only at very high density, 
which, as we will see, is actually not reached in the corre¬ 
sponding NS structure. 

Once the EOS of the nuclear matter which is present 
throughout the NS is known, one can use the Tolman- 
Oppenheimer-Volkoff llobl - l^ equations for spherically sym¬ 
metric NS: 


dp 

dr 

dm 

dr 


Gme (1 + p/e)(l +47tr^p/m) 
r^ 1 —IGm/r 

= 4;rr^e, 


( 21 ) 

( 22 ) 


where G is the gravitational constant and m{r) is the enclosed 
mass within a radius r. Given a starting density e^, one inte¬ 
grates these equations until the surface r = R, and the grav¬ 
itational mass is obtained by Mq = m{R). The EOS needed 
to solve the TOV equations is taken from the nuclear mat¬ 
ter calculations as discussed above for the liquid-core region 
and match ed w ith the crust EOS, which has been taken from 
Refs. lIMMl- This matching occurs at about two thirds 
of the saturation density, where the EOS analytical fits of 
Eq. (fPl) are still accurate. 

As is well known, the mass of the NS has a maximum value 
as a function of radius (or central density), above which the 
star is unstable against collapse to a black hole. The value of 
the maximum mass depends on the nuclear EOS, so that the 
observation of a mass higher than the maximum mass allowed 
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FIG. 11. (Color online) Neutron star mass as a function of radius (left 
panel) or central baryon density (right panel) for several EOS. Thin 
lines indicate results obtained with purely nucleonic EOS, whereas 
thick lines show results for several EOS including hyperons. See text 
for details. 


by a given EOS simply rules out that EOS. The fss2 EOS 
gives slightly different maximum masses for the gap choice 
and continuous cases, in line with their different stiffness at 
high density. This gives a range of uncertainty for the maxi¬ 
mum mass, which encompasses the large st mass observed up 
to now, which is (2.01 ±O.O4)M0 llO^ . This is illustrated 
in Eig. [TTl where the relation between mass and radius (left 
panel) or central density (right panel) is displayed in compar¬ 
ison with the other considered EOS (thin lines). The observed 
trend of the central density for all the EOS is consistent with 
the corresponding P{p) relation displayed in the central panel 
of Eig. [To] As expected, with increasing incompressibility, the 
NS central density decreases for a given mass. 

In the end we illustrate the so-called “hyperon puzzle” with 
the fss2 model. Eig. [TT] (thick lines) shows the effect of allow¬ 
ing the appearance of hyper ons in beta-stable matter within 
our BHF approach lll03lll04|] . Two different nuc leon- hyperon 
(NY) interactions, th e Nij megen NSC89 model lilOSll and the 
recent ESC08 model lllObll are considered, and combined with 
the fss2 NN potenti al in the approximate way explained in 
more detail in Ref. lll03l] . namely the purely nucleonic BHF 
energy density obtained with the fss2 is combined with the 
hyperonic contribution to the energy density evaluated with 
either the NSC89 or the ESC08 interaction, but together with 
the Argonne Tig potential plus nucleonic TBF. In this way the 
intermediate states in the NY Bethe-Goldstone equation are 
treated approximately, but t he ov erall error of the global re¬ 
sults is expected to be small lll03l] . 

Fig. [TT] demonstrates that under these assumptions the NS 
maximum mass is practically insensitive to the choice of the 
NN interaction, but determined by the NY interaction. This 
is due to a well-known compensation mechanism that can be 
clearly seen in the left panel of the figure: The slightly stiffer 


fss2 gap-choice model (dashed black and blue curves) causes 
an earlier onset of hyperons and a stronger softening than the 
fss2 continuous-choice model (solid curves). In any situation 
the maximum mass is much smaller than current observational 
values. 


V. SUMMARY AND DISCUSSION 

We have derived the nuclear matter EOS within the BBG 
approach up to the THE level of approximation, starting from 
an NN interaction based on quark-quark and quark-meson in¬ 
teractions. An intrinsic uncertainty in the approach is related 
to the choice of the auxiliary s.p. potential. Within this un¬ 
certainty the saturation point is well reproduced without any 
additional parameters with respect to the interaction, that is 
able also to reproduce the binding of three and four nucleon 
systems. At higher energy the interaction should be improved 
in some channels, in particular the ^P 2 -^F 2 , which is relevant 
for the high-density part of the EOS, and therefore for NS. 

The symmetry energy as a function of density up to satura¬ 
tion, its value and slope at saturation, and the incompressibil¬ 
ity of symmetric matter at saturation compare favorably well 
with the phenomenological constraints. Above saturation the 
EOS is compatible with the flow data in HIC at intermediate 
energy, up to about four times saturation density. 

As already discussed, a warning about these calculations is 
the observation that two and three hole-line diagram contribu¬ 
tions become comparable at higher density, which puts some 
doubt on the convergence of the BBG expansion. However, 
this is mainly due to the behavior of the two hole-line con¬ 
tribution, which saturates or even decreases at higher density, 
due to the compensation between positive and negative con¬ 
tributions. This indicates that the degree of convergence can¬ 
not be estimated in a straightforward way. Within the present 
many-body treatment this is probably the main source of un¬ 
certainty in the results. 

However, let us notice that up to few times saturation den¬ 
sity the EOS calculated with the continuous and gap choice 
agree very well. In particular for symmetric matter the agree¬ 
ment extends up to the maximal density used in NS calcula¬ 
tions. This fact can be considered a good indication for the 
convergence of the expansion, because this agreement would 
be exactly true if convergence is indeed reached. For similar 
reasons the mass-radius relationship in NS are quite close for 
the two choices, which gives support to the validity of these 
results. 

The EOS can be considered relatively soft, but despite that 
the NS maximum mass is compatible with the curre nt ob - 
served NS maximum mass of about 2 solar masses ll02|] . 
Phenomenology seems then to validate this microscopically 
derived EOS, at least up to few times saturation density. 

However, there are some theoretical caveats to be consid¬ 
ered. It can be expected that quark matter appears in the cen¬ 
ter of massive NS. To describe these “hybrid” NS one needs to 
know the quark matter EOS. It turns out that many models for 
the deconfined quark matter produce a too soft E OS to sup - 
port a NS of mass compatible with observations IIOTITI 13n . 












13 


The quark-quark interaction in the deconfined phase must be 
then repulsive enough to stiffen the EOS, and indeed, with a 
suitable quark-quark interaction, mixed quark-nucleon matter 
can have an EOS compatible with two solar masses or more 

iniini. 

An additional problem arises if strange matter is introduced 
in the NS matter. It turns out that BHE calculations using 
realistic hyperon-nucleon interactions known in the literature 
produce a too soft NS matter EOS and the maxi mum mas s 
is reduced to well below the observational limit iflolflol . 
Although the hyperon-hyperon interactions and in particular 
hyperonic TBE are poorly known, these results pose a “hy- 
peron puzzle”. EOS b ased on re lativistic mean field models 
can solve the problem ill IhUl isll with a proper choice of pa¬ 
rameters. Also modifications of the hyperon-nucleon inter¬ 
action, including thr ee-body fo rces, could provide a remedy 
for the too soft EOS ill 19[ll2Cill . All these methods introduce 
quark-quark or nucleon-hyperon and hyperon-hyperon inter¬ 
actions that stiffen enough the EOS at high density. 

In this respect it would be of particular relevance to per¬ 
form BBG calculations up to THE level with the quark-model 
baryon-baryon interaction fss 2 extended to the strange sector 
ll44ll . This difficult problem must be left to a future long-term 
project. 
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Appendix A: Deuteron wave function 

In this appendix, we basically follow the notation of 
Ref. lH^. Eirst we solye the Lippmann-Schwinger ROM 
equation for the deuteron 

/ dqq^Vee,{k,q)fe,{q). 

(2ny yjo 

(Al) 

Here we use the nonrelatiyistic expression = J^/Mn as the 
deuteron binding energy. The relatiyistic correction is of order 
Ej/Mn, which corresponds to a few keV difference if^ . The 
total waye function is 

'I'f (fc) = [Mk)%\'^{k)+f2{k)^2\^{k)] Co° , (A2) 

where '3^l^{k) are spin-spherical harmonics and the isospin 
function is denoted by . In coordinate space, we haye 

= [(MoWA)^r(^)+(M2WA)^r(^)] Co ■ 

(A3) 


They are related by the Eourier transform 

ut{r) = \/Tf% y f dkiyrji{kr)fi{k) . (A4) 

Jo 


TABLE VIII. The range parameters Yj and coefficients Cj and Dj for 
the parametrized deuteron wave function. The values in parentheses 
are calculated from the boundary conditions. 


J 

Yi 

Cj 

Di 

1 

0.2314 

0.8803 

0.2218 X 10^* 

2 

1.1314 

-0.2235 

-0.4637 

3 

2.0314 

-0.2501 

0.1776 

4 

2.9314 

-0.1658 xlO^ 

-0.1294 X 10^ 

5 

3.8314 

0.1694 X 10^ 

0.1192X 10^ 

6 

4.7314 

-0.9160 X 10^ 

-0.5884 X 10^ 

7 

5.6314 

0.2592 X 10"^ 

0.1611 X 10"^ 

8 

6.5314 

-0.4095 X 10'^ 

-0.2503 X lO'* 

9 

7.4314 

0.3678 X lO'* 

(0.2214 X lO'*) 

10 

8.3314 

-0.1766 X 10^ 

(-0.1043 xlO^) 

11 

9.2314 

(0.3530 X 10^) 

(0.2040 X 10^) 


The normalization is 

Y [ druj{r)= Y [ dkk^fj{k) = \. (A5) 

1=0,2-'^ f=0.2-'0 

The deuteron wave functions are parametrized in the fol¬ 
lowing way as in Refs, ll^ and if^ 

huik) = E { % } ViAi? < A • 'A® 

n r Cjc-yi'' 

uia{r) = Y\ + ^ + 

(A7) 

The range parameters are chosen as 7 ; = 7 + {j — l)7o with 
7 o = 0.9 fm^* and « = 11. 

Eor r the deuteron wavefunctions have the form uo{r) 
-t-and U 2 {r) -i [l-f 3 /( 7 ir)-f 3 /( 7 ir) 2 ], 

where A^ = Ci and Ao = D\. The asymptotic D-state to S- 
state ratio is given by 77 = As lAjj. The boundary conditions 
Mo(r) r and U 2 {r) —as r —> 0 le ad to one constraint for 
the Cj and three constraints for the Dj CH. These constraint 
for the last Cj and the last three Dj are explicitly written in 
Eqs. (C.7) and (C. 8 ) of Ref. ll^ . The values of jj, Cj, and Dj 
are listed in Tab. IVIIII The accuracy of the parametrization is 
characterized by 

|^“flfkA[/o(k)-/oa(k)f| ^ =1.9 x10-A (A 8 ) 

dkk^[f2{k)-f2a{k)fy =2.3x10-A (A9) 

The quadrupole moment Q^, the root mean square radius 
Rd, and the D-state probability Pj) are calculated using 
Eqs. (C.16)-(C.18)of Ref. if^ . respectively. 


r £ = 0 
{£ = 2 ■ 
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